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ABSTRACT 

The spectral energy distribution analysis of very young unresolved star clusters challenges our 
understanding of the cluster formation process. Studies of resolved massive clusters in the Milky 
Way and in the nearby Magellanic Clouds show us that the contribution from photoionized gas is 
very important during the first Myr of cluster evolution. We present our models which include both a 
self-consistent treatment of the photoionized gas and the stellar continuum and quantify the impact of 
such nebular component on the total flux of young unresolved star clusters. A comparison with other 
available models is considered. The very young star clusters in the SBS 0335-052E dwarf starburst 
galaxy are used as a test for our models. Due to the low metallicity of the galactic medium our models 
predict a longer lasted nebular phase which contributes between 10-40 % of the total near infrared 
(NIR) fluxes at around 10 Myr. We propose thus a possible solution for the observed flux excess in the 
6 bright super star clusters of SBS 0335-052E. Reines et al. showed that the observed cluster fluxes, 
in the red-optical and NIR range, sit irreconcilably above the provided stellar continuum models. We 
find that in the age range estimated from the Ha emission we can explain the red excess in all the 6 
super star clusters as due to nebular emission, which at cluster ages around 10 Myr still affects the 
NIR wavebands substantially. 

Subject headings: galaxies: individual (SBS 0335-052E) - galaxies: dwarf - galaxies: starburst - galax- 
ies: star clusters 



1. INTRODUCTION 

SBS 0335-052E, at a distance of 54 Mpc (according 
to NED), is one of the least chemically evolved galaxies 
in the local universe (Z=l/40 Z©, Izotov et al. 1990; 
Izotov et al. 1997; Kunth & Osthn 2000; Papaderos et 
al. 2006). In the 1990s, the first Hubble Space Telescope 
(HST) images of this irregular dwarf galaxy (Thuan et 
al. 1997) revealed that the main starburst region was 
dominated by 6 young super star clusters (SSCs). In 
Figure 1 we show the position of the clusters as identified 
by Thuan et al. (1997) using the nomenclature of Reines 
et al. (2008b, hereafter R08). The clusters are located 
around a supernova cavity which has probably caused 
their formation, with star formation propagating from 
north to south (Izotov et al. 1997). 

Thompson et al. (2009, hereafter T09) observed bright 
Paa emission in the galaxy and identified three new and 
actively star forming regions, S3-Paa, S7 and S8 (in- 
dicated in the continuum subtracted Ha image in Fig- 
ure 1). These new areas were presented as places of less 
clustered star formation, where young stars are ionizing 
the gas in situ. Similar star forming regions were also 
observed around SSCl and SSC2. 

Vanzi et al. (2000), found that a hot dust component is 
needed to explain the red H — K color observed globally 
in the galaxy. However, the very low visual extinction 
across the SSCs, which show prominent Br7 emission, 
implies that the clusters are quite dust free. Interestingly, 
Hunt et al. (2001), using imaging and spectroscopy data 
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at 2 — 4 /im, found evidence of a deeply embedded cluster 
with a visual extinction Av > 15 mag in the line of sight 
of the two young SSCl and 2. They excluded that the IR 
fiux was produced by the two optical bright SSCs because 
of the high discrepancy between the extinction derived 
from the IR and from the optical data {Ay ^ 0.55 from 
the optical line ratios). 

A different scenario was presented by ROB and John- 
son et al. (2009) to reconcile the apparently discrepant 
extinction estimates in the knot containing SSCl and 2. 
ROB inferred ^ 40% of Lyman continuum photon leakage 
from the two clusters which is ionizing the extended Ha- 
bright region around the knot. This could be possible 
if the medium surrounding the two clusters is clumpy. 
In this case the optical light would come through the 
"holes" suffering only low extinction, while in the NIR 
would transmit also the signal produced by the dense 
dust clumps, heavily extinguished. 

Using Ha equivalent width EW(Ha) measurements, 
ROB deriveded cluster ages between 3 and 15 Myr. The 
6 SSCs also displayed a flux excess at wavelengths larger 
than BOOO A, impossible to reconcile with the stellar 
continuum models used to fit the cluster spectral en- 
ergy distribution (SED). ROB explored several possibil- 
ities to explain the origin of the red excess. The excess 
at ~ O.B//m was explained as the dust photoluminescence 
phenomenon commonly referred to as extended red emis- 
sion (ERE, see Witt & Vijh 2004). The ERE has been 
observed in star forming regions like 30 Doradus (Darbon 
et al. 199B) and in Hii regions of local galaxies like NGC 
4826 (Pierini et al. 2002), and is characterized by an ex- 
tended emission feature between 7000 and 9000 A. ROB 
attributed the fiux excess at wavebands between 1.6 and 
2.1 /zm to two different mechanisms. For the youngest 
clusters SSCl and SSC2, hot dust (- BOO K) emission 
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was advocated to explain the rise of the IR continuum. 
The remaining older clusters showed infrared colors con- 
sistent with those of red super giants (RSGs). In fact, 
if the models do not contain a reliable treatment of the 
RSGs (Langon et al. 2008 and references therein) this 
could produce an apparent NIR excess in the observed 
SEDs at cluster ages around 10 Myr. 

The clusters are very young (3 — 15 Myr). This age 
range is quite complex, due to the rapid dynamical evo- 
lution that the clusters experience (see Portegics Zwart 
et al. 2010 for a review). At such young stages, the 
contribution from photoionized gas in the broadband in- 
tegrated fluxes can be substantial ( Bergvall & Ostlin 
2002; Anders & Fritze-Alvensleben 2003; Zackrisson et 
al. 2008; Reines et al. 2010). In the starburst galaxy 
Haro 11, we found a peak of cluster formation only 3.5 
Myr old (Adamo et al. 2010). The analysis of such young 
clusters requires models that include both the stellar con- 
tinuum and photoionized gas emission contributions (see 
Zackrisson et al. 2001, and Adamo et al. 2010). How- 
ever, we also found that a subsample of clusters in Haro 
11 dysplays a flux excess above models including pho- 
toionized gas emission at wavelengths longward of 8000 
A. The observed discrepancies could be explained if other 
mechanisms, not included in the models and related to 
the young ages of the clusters, are contributing at red- 
der wavebands. Hot dust, a considerable fraction of prc- 
main sequence stars (PMSs) and young stellar objects 
(YSOs), are possible mechanisms causing the red excess 
in Haroll. Previous cases of NIR excess observed in 
still embedded star clusters of dwarf starburst galaxies 
was explained by several authors as caused by hot dust 
(Cresci et al. 2005; Cabanac et al. 2005). Observations of 
very young star clusters in the Small Magellanic Clouds 
have revealed complex environments composed by bright 
main sequence stars, but also pre-main sequence objects 
(PMS) and young stellar objects (YSOs) only observable 
in the NIR (Carlson et al. 2007; Gouliermis et al. 2010). 
A significant fraction of YSOs was suggested as the ori- 
gin of the red excess in the starburst galaxy NGC 253 
(Fernandez-Ontiveros et al. 2009). 

The NIR excess in young star chistcrs has been ob- 
served in several nearby starburst galaxies. Here, we 
explore to what extent the red excess reported in SBS 
0335-052E by R08 can be explained by nebular emission, 
or other mechanisms are needed. 

In Sections 2 and 3 we present the photometric data 
and our synthetic evolutionary models which we use to 
perform new fits to the cluster SEDs as described in Sec- 
tion 4. In this section we also test the validity of our 
fits and attempt an estimation of the sizes of the emit- 
ting Hii regions. In Section 5 we compare our results 
with previous studies of the galaxy. Our conclusions are 
summarized in Section 6. 

2. DATA SAMPLE 

The data reduction and photometry were carried out 
and described in detail in Section 3 of R08. In the present 

work we will use their photometry (published in their 
Table 1) to test our models and investigate the origin of 
the red excess in SBS 0335-052E. 

In Figure 1 we show two of the science frames of SBS 
0335-052E, previously published by Ostlin et al. (2009). 
We illustrate the position of the clusters as in R08 and 



the star forming regions identified by T09. The photo- 
metric annuli around the clusters have a radius of 0.15" 
(~40 pc at the adopted distance). This corresponds to 
the size set by R08 to perform aperture photometry. The 
six SSCs are bright and visible in all the UV, Optical and 
IR i?5r frames. The continuum frame in Figure 1 shows 
that the positions of the compact clusters are clearly de- 
fined and that the used aperture size is suitable to avoid 
contamination from the neighboring diffuse regions. On 
the other hand, in the right frame, we see that Ha emis- 
sion is quite strong at the location of the SSCl and SSC2 
which makes it difficult to distinguish between the com- 
pact and the less clustered surrounding regions. The re- 
maining 4 star clusters do not show Ha fluxes as strong 
as the ones produced in the nearby very young active 
star forming regions. 

2.1. New estimation of the EW(Ha) in the clusters. 

The use of local sky annuli, like the ones used by R08, 
could oversubtract flux from the sources, producing an 
underestimation of the final EW(Ha). Due to the com- 
plexity of the environment around the clusters, we de- 
cided to re-estimate the EW(Ha) across the clusters us- 
ing the same aperture radius as R08, but a constant value 
for the background as explained below. 

The reduction of science frames shown in Figure 1 
and software arc described in more detail in Hayes et 
al. (2009), here, we briefly summarize the procedure. 
The continuum subtraction of the Ha frame {FR656N) 
is non-standard, and utilizes custom software. This soft- 
ware performs spatially resolved, pixel-by-pixel SED flt- 
ting using five HST broadband filters between the FUV 
and optical domain, all of which were expressly chosen to 
be free from contamination by nebular lines. From the 
best-fitting synthetic stellar population, it computes the 
expected flux due to continuum processes in the FR656N 
narrowband filter. This is then subtracted from the ob- 
servation itself, and maps of the pure 6563 A continuum, 
and continuum-subtracted Ha are output. 

Photometry was carried out in both frames, at the 
same position, using the same aperture radius (0.15") as 
in R08. Because the Ha emission has a complex and ex- 
tended morphology, we preferred to subtract a constant 
mean value of the sky background instead of a local sky 
annulus located around the position of the cluster, in 
both continuum and emission frame. The mean sky val- 
ues were estimated in both frames, in the same position, 
and in a region that appeared free of Ha emission and 
continuum flux. No aperture corrections to the deter- 
mined fluxes were applied since there are no point-like 
sources in the reduced frames and simulated PSFs would 
introduce additional uncertainties. For this reason, it is 
possible that the new estimated EWs could be a lower 
limit to the real values. However, the EW is a ratio, 
so if line and continuum have the same distribution (an 
assumption that is valid for high EWs) the aperture cor- 
rection would be cancelled. 

The EWs were estimated by the ratio between the Ha 
emission flux multiplied by the filter width and the flux 
in the continuum. In Table 1, measured Ha flux and 
EWs are listed for all the clusters. In general, the ob- 
served EW(Ha) found by R08 are smaller than the ones 
estimated here. The main differences between the two 
analyses reside in the reduction of the frames, and in 
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Fig. 1. — Ho continuum and emission frames of SBS 0335-052 (Ostlin et al. 2009). The position of the compact SSCs is showed by the 
red circles. The radius of the circles is 0.15". The new star forming regions with bright Ha emission are indicated in the right frame using 
the nomenclature by T09. 



the performed sky background subtraction. The contin- 
uum frame, in our analysis, is better constrained than 
the one produced by a simple interpolation between the 
two nearest frames to Ha. 

3. DESCRIPTION OF THE SPECTRAL EVOLUTIONARY 
MODELS 

To produce self-consistent models which include the 
contributions from both stars and photoionized gas, we 
use the stellar population SEDs predicted by Starburst99 
(Leitherer et al. 1999) as input to the Cloudy photion- 
ization code (Ferland et al. 1998). The procedure closely 
follows that described by Zackrisson et al. (2001), in 
which Cloudy is called for every stellar population time 
step to ensure that the temporal evolution of the neb- 
ulae is treated self consistently. This technique results 
in combined SEDs that are more realistic than the ones 
generated when using the nebular components produced 
by Starburst99 itself. The latter contain only nebular 
continuum (produced by free-free and free-bound emis- 
sion), whereas the models we use include both nebular 
continuum and emission lines (produced by bound-bound 
mechanisms). In our models, the SB99 stellar population 
SEDs are produced using the Padova tracks (Bertelli et 
al. 1994) and assume a metallicity of Z = 0.0004, an in- 
stantaneous burst of star formation (resulting in a single- 
age population) and a stellar initial mass function with 
a double power-law (dN/dM cx Af~", with a = 1.3 for 
M = 0.1-0.5 Mq and a = 2.3 for M = 0.5-100 Mq). 
The nebular SED assumes a spherical, constant-density 
and ionization-bounded nebula with a gaseous metallic- 
ity identical to that of the stars. The gas filling factor 
/ (which quantifies the porosity of the nebula), the hy- 
drogen number density nn, and the gas covering factor c 
(which quantifies the amount of Lyman continuum leak- 
age due to dumpiness in the medium or presence of dust) 
are all explored within reasonable boundaries. We gen- 
erated a grid of models, with different combinations of 
values for / = [0.01,0.001], uh = [10^10^] cm'^, and 



c = [0.1,0.2,0.3,0.4,0.5,1.0] (c = 1 means that all the 
Lyman continuum photons ionize the gas; c = 0.1 that 
only the 10 % of them photoionizes the gas, i.e. that 90 
% of photons arc leaking). 

In the Figure 2, we show how the integrated spectrum 
of a 10^ Mq cluster changes as a function of the age dur- 
ing the first 15 Myr, using our models with / = 0.01, 
TiH = 10^ cm~^, and c = 1. For comparison we plot the 
corresponding Starburst99 stellar continuum spectrum 
(red dotted line) and the model with nebular continuum 
included (blue line). In the very early stages, around 
3 Myr, the contribution from the photoionized gas is so 
strong that the whole spectral range from the near-UV to 
IR is affected. The stellar continuum is significantly be- 
low. The continua of the two models including the neb- 
ular contribution are in good agreement. However our 
models, implemented with Cloudy, include also emission 
lines, very strong at this age. At 6 Myr, we observe that 
the stellar continuum and the spectrum of the model in- 
cluding both stellar and nebular components are similar 
in the near-UV and blue-optical range. The contribution 
from the line emission, however, is very important in the 
Balmer region. From the red-optical to the near-IR range 
the difference in fluxes between models with and without 
nebular treatment is still substantial. Around 10 Myr, 
fluxes at wavelengths longward 1.0 /im still receive a non 
negligible fraction of flux from the gas. This is evident 
using both SB99 (stellar and nebular continuum) and 
our models. At shorter wavelengths, only the integrated 
fluxes across the emission lines are affected. The nebular 
component is almost gone at around 15 Myr, with very 
small differences between the 3 spectra. 

In Figure 3 we compare the fraction of nebular emis- 
sion contributing to the total flux as estimated by our 
models and by the SB99 ones. We quantify these frac- 
tions as a function of the cluster age for the HST fllters 
used in the analysis of SBS 0335-052E (R08). As already 
showed in the previous flgure, the only significant differ- 
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Fig. 2. — Evolution of the spectral model for a IO^Mq cluster. In black we show the model including a full treatment on the nebular 
emission and continuum, in red the stellar continuum model, and in blue the stellar and nebular continuum only (the last two models are 
outputs from Starburst99). The fluxes are in arbitrary units. The wave range is between 0.1 and 2.3 fim. 



ence between the two models is the presence/ absence of 
emission Unes. On the top panels, we clearly see the drop 
in the contribution of the photoionized gas already at 
4-6 Myr in the NUV (.4C5'/F220W), U (^C5/F330W), 
B (^C5/F435W), and V {ACS/F550M) filters. During 
this age range, the difference between the two models 
is of < 5 % and disappears at older ages, if the filter 
does not transmit important lines (i.e. NUV, U, and V) . 
On the other hand, the B filter shows in our models a 
significantly higher fraction of photoionized gas contri- 
bution, due to the presence of the Balmer lines. Around 
10 Myr the B band still transmits 10 % of the nebular 
flux which is mainly contained in the lines (see Figure 2 
at ~ 0.4 — 0.45 /im). In the remaining UV and opti- 
cal bands, only 5 % of the integrated flux is produced 
by the nebula. In the bottom panels we show the NIR 
fllters. Again the discrepancy between our models and 
the SB99 ones is higher in filters which transmit strong 
lines (mainly in K). As already noticed before in the NIR 
range, at 10 Myr, the nebula still supplies 10 % of the to- 
tal flux in the optical I ( WFPC2/F79m) band, 20 % in 
the H (7V/C^/F160W) band and in the Paa-continuum 
narrow fiher (Af/C^/F187N), and up to 40 % in the K 
(7V/C^/F205W) band. We also notice that at a given age 
the fraction of flux produced by the ionized gas increases 
as function of the wavelength, and shows a slower decline 
towards older ages. 

As further test, in Figure 4 we plot the evolutionary 
tracks in different combination of IR and optical colors. 
As already observed in Zackrisson et al. (2001), during 
the first 10 Myr the predicted colors by models which 
include contribution from photoionized gas are quite dif- 
ferent from the ones with only stellar continuum. We 



also plot, in the color diagram, the position of the six 
SSCs of SBS 0335-052E. We clearly see that the cluster 
IR colors are much redder than the predictions provided 
by the stellar continuum only or the stellar and nebular 
continuum tracks (SB99), but in fairly good agreement 
with our models. The position of the clusters in the 
bottom panel is shifted toward red F550M-F791W color 
{V — I > 0). This displacement is due to a small flux 
excess in I band, which is around 0.2 — 0.4 mag above 
the best fitting models. This issue will be discussed in 
the following sections. 

4. CONSTRAINING THE ORIGIN OF THE RED EXCESS 
4.1. A new cluster SED fit from the NUV to the NIR 

The SED fit made by R08 was limited to the near- 
UV and blue-optical fihers (y4C5/F220W, ^C5'/F330W, 
^C5/F435W, ACS/Ybbm). The Starburst99 spectral 
evolutionary models, including only the stellar contin- 
uum component were used together with the LMC ex- 
tinction curve (Misselt et al. 1999, and references therein) 
to estimate ages, masses and extinctions of the clusters. 

In Section 3, we argued that the photoionized gas is a 
crucial component for the models used to estimate the 
physical properties of the young clusters. In this section 
we explore the possibility to explain the red excess ob- 
served by R08 using these new models. The young clus- 
ters in SBS 0335-052E are indeed in the the age range 
where the contribution from the photoionized gas sur- 
rounding the clusters is expected to be important. The 
amount of fiux produced by the gas and contributing to 
each of the cluster integrated fiuxes is a function of the 
total amount of ionizing UV-flux produced by the mas- 
sive stars in the cluster. The EW(IIa) is a direct mea- 
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Fig. 3. — Relative contribution of the nebular emission as func- 
tion of the cluster age to the total flux transmitted in the HST 
filters used in this analysis. The SB99 models including stellar 
and nebular continuum are showed with a dashed dotted line. Our 
models are drawn by solid lines. The filters are indicated in the 
plots. A description of the models is given in section 3. 

surement of the strength of the ionizing flux and very 
sensitive to the rapid evolution of the massive stars dur- 
ing the first 15 Myr of the hfe of the cluster. R08 com- 
pared the measurement of the EW(Ha) with predictions 
made by SB99 to derive the ages of the clusters. We per- 
formed a similar estimation using the new measurements 
of EW(Hq!) of each SSC and the predictions provided by 
Cloudy in our models. The ages constrained by R08 and 
ours are reported in Table 1 (fourth and fifth column, re- 
spectively). The ages we derived are in good agreement 
with what already found by R08 and by Johnson et al. 
(2009) for the two bright radio sources SSCl and 2. 

Keeping fixed the age, we produced a fit of the ob- 
served integrated fiuxes of each cluster from UV to IR 
bands with the fluxes provided by our models at the cor- 
responding ages. Our fit algorithm is described in detail 
in Adamo et al. (2010). Internal extinction was treated 
as a free parameter, and we used the Calzetti attenua- 
tion law (Calzetti et al. 2000), allowing extinction vary- 
ing from E{B - V) = 0.0 to 3.0 with a step of 0.01. In 
the previous fit, R08 used the LMC extinction law. We 
showed however, in Adamo et al. (2010), that there are 
not substantial deviations in the estimated ages, masses, 
and number of clusters affected by a red excess, when a 
Calzetti or a LMC law are used. The mass of the clus- 
ters was derived from the normalization between best- 
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Fig. 4.— The FlS7N-F20bW color versus the F330VF-F435VK 
and the redder i^550M — ones in two color diagrams. The 

evolutionary track provided by our model is represented by a solid 
line. We included Starburst99 evolutionary tracks with only stellar 
continuum (red thick dotted line) and stellar plus nebular contin- 
uum (blue thick dashed line) as well. Positions in the tracks corre- 
sponding to different ages are labelled. The colors of the six SSCs 
in SBS 0335-052E are overplotted. The arrow shows in which direc- 
tion the correction for an visual extinction of = 1.0 mag moves 
the position of the clusters in the two diagrams. Photometric error 
bars are included. 

fitted model and observed data. The models that pro- 
duce smaller residuals are the ones with: / = 0.01,c=l, 
and uh — 10^ cm^^. The fits are shown in Figure 5, and 
the masses and A\/ are given in the Table 1. We noticed, 
however, that the models with same / and c values but 
lower hydrogen densities (n// = 10^ cm~'^) produce fits 
with residuals only slightly bigger than the best models. 
The constrained extinction range is in good agreement 
with the values found using optical and NIR line ratio 
(e.g. Vanzi et al. 2000). The values are also in agree- 
ment with the recovered extinction ranges in other star- 
burst environments like the Antennae system (Mengel et 
al. 2005) and M51(Bastian et al. 2005). The estimated 
cluster masses for the two youngest clusters, SSCl and 
2, are significantly lower than the ones determined by 
R08. They are instead, in agreement with the estimates 
by Johnson et al. (2009) derived from the thermal radio 
emission of the two young SSCs, under the assumption 
of no leakage of ionizing photons from the nebula sur- 
rounding the young cluster. The remaining clusters have 
estimated masses in agreement with the values previously 
obtained by R08. During the first few Myr, where the 
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photoionized gas emission is dominant, we clearly see 
that the cluster mass can be overestimated by a factor 
of 2 — 3, if the used models do not include nebular con- 
tribution. 

We also repeated the least-square fit done by R08, ex- 
cluding the I band and IR data ( WFPC2l¥l%m, and 
NIC2 F160W, F187N, and F205W), and using models 
with only stellar continua. In this case, age was treated 
as a free parameter. The fits with the smallest residuals 
are shown in Figure 5 with red data points and red spec- 
tra. As already found by R08, the fit to the blue optical 
photometry using only stellar continua produced a quite 
evident displacement at wavelength longward of 8000 A. 
Nevertheless, when nebular continuum and line emission 
are included in the models, we obtained a good fit for 
all the clusters and in all the passbands. Inside the pho- 
tometric errors, the best fitting model (black squares) is 
able to reproduce the SED shapes quite well. We noticed, 
however, a persisting small excess (between 0.2 and 0.5 
mag above the models) at 8000 A (I band) which was ob- 
served in the bottom panel of Figure 4 . We will discuss 
this feature in Section 5. 

We produced, finally, a third set of fits from UV to 
IR, using our models, where we let the age be a free 
parameter. The two youngest SSCl and SSC2 had a 
best fitting age of 3 Myr in perfect agreement with the 
EW(Ha) determination. The SSC3, SSC4, SSC5, and 
SSC6 had best fitting ages lower than the values deter- 
mined using EW(Hq!): 6, 9, 9, and 7 Myr, respectively. 
However, investigating the 68.3 % confidence levels, we 
constrained boundaries to the best fitting ages (Table 2). 
Wc observe that inside the imccrtaintics produced by a 
free SED fit the solutions recovered using the observed 
EW(Ha) are in good agreement in all the cases except 
in SSC6 (marginally agreement). The estimated masses 
and extinction didn't change significantly (see Table 2). 
Therefore wc do not discard any of the two results and 
consider the ages determined with the two methods as a 
valid range around the real age of the clusters. 

4.2. Testing the validity of the models 

Although our models can reproduce the shapes of the 
cluster SEDs, we need to test the consistency between the 
provided models and the physical conditions of the clus- 
ters. The models include the full amount of emission flux 
coming from the photoionized gas. This assumption im- 
plies that all the ionized gas (e.g. the Hll region) around 
the cluster is spatially contained in the used photometric 
aperture size. Using the nebular properties produced by 
the best fitting model, and cluster physical parameters 
published in previous works, we produced estimates of 
the radii of the Hii regions around the clusters, as ex- 
pected from the model. The radii were computed using 
the formula: 

under the assumption of a pure hydrogen cloud. 

The Q Lyc is the rate of ionizing Lyman continuum pho- 
tons produced by very young stars and changes as a func- 
tion of cluster age. If the medium around the stars is 
clumpy or dusty, it is possible that a fraction of the ion- 
izing photons do not ionized the gas. R08 discussed the 



possibility of a clumpy medium around the two youngest 
SSCl and 2, and estimated a leakage of ionizing photons 
of ~ 40%. Our best fitting model had a covering factor 
c = 1, implying that all the Lyman continuum photons, 
produced by the young massive stars, ionize the gas. To 
verify this result we compared the provided Q Lyc by the 
Starburst99 model with the values of QLyc derived from 
the measured Ha fluxes (extinction corrected). The ob- 
served Qhyc was estimated using the Relation (1) in R08, 
which produces lower limits to the real values. The f(Ha) 
listed in our Table 1, were first corrected for the corre- 
sponding extinction values, and then converted into H/3 
luminosities. The resulting Qhyc were scaled for the the 
mass of the model (lO^M©). We show the calculated 
values of the QLyc for the 6 SSCs in Figure 6. For the 
two youngest clusters the estimated QLyc values are in 
agreement with the model predictions, and do not show 
any significant leakage of ionizing photons. We observe 
in the older clusters a possible escaping fraction of 20—30 
%, consistent with the more evolved phases the clusters 
are experiencing. In general, we consider the assumption 
of c = 1 a reasonable approximation. 

To estimate the radii of the Hii regions of each cluster 
we used the QLyc predicted by the Starburst99 stellar 
continua, scaled for the mass of the cluster. In Table 1 
we also show (into brackets) the values of the radii if the 
estimated QLyc from f(HQ;) woifld be used instead. 

Studies of massive star forming regions in our galaxy 
(Churchwell 2004, Conti & Blum 2002, among many oth- 
ers), have shown qiiitc complex environments with den- 
sity gradient going from nn ~ lO'^ cm~^ in pre-stellar 
cores to diffuse surrounding regions with uh ~ 10^ cm^'"', 
and ultra-compact Hii regions {uh ~ 10^ cm^'*), all co- 
existing in the same system. Thompson et al. (2006) 
presented observational evidence of very compact Hll re- 
gions in the position of the two yoimg SSCl and SSC2. 
On the other hand, using two-dimensional spectroscopy, 
Izotov et al. (2006) found values of ur across the SSCs of 
a few hundreds and EW(Ha) below 1500 A even in the 
brightest Ha regions of the galaxy (SSCl-l-2). We no- 
ticed, however, that the spectra are produced in regions 
approximate 10 times bigger (0.52" xO.52" ) than the ones 
used here (radius of ~0.15"), causing a smoothing of the 
values over low density and diffuse regions around the 
clusters. Recently using continuum radio data, Johnson 
et al. (2009) estimated for SSCl and 2 a n// ~ 4-7 x 10^ 
cm^"^ in agreement with the value of Hh = 10* cm~^ we 
derived from the fitting model. 

In order to estimate the radii of the Hll regions for 
the two youngest clusters we used hr = 10* cm~^. In 
the previous section, the model with lower uh produced 
also a good fit to the data. By keeping all the other 
parameters fixed, wc checked whether a lower value of 
Hh = 100 cm^'^ produced very different fiux values in 
the optical I band and in the NIR wavebands for ages 
between 6 and 15 Myr. We found that the difference 
between the two models was less then 1%. Finally, we 
decided to use nn = 10^ cm~'^ for SSC3, which has an 
intermediate age, and a density of few hundreds, for the 
older three SSCs. 

The filling factor, f, was fixed to 0.01. We set the re- 
combination coefficient, a^, for Case B recombination 
(Osterbrock & Ferland 2006), to 1.43x10"" at an elec- 
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TABLE 1 

Properties of the 6 SSCs. SED at fixed age for mass determination. Age is determined 

TTiRorcii EWVIIn V 



Ids 


f(Ha) 


EW(Ha) 


Age (R08) 


Age'' 


Mass 


Av 


f-Hii 




[erg s"-*- cm~^] 


A 


Myr 


Myr 


Mq 


mag 


pc 


SSCl 


5.14x10-1* (1%) 


3100.0(134.4) 


< 3.3 


3.0 


4.7x105 


0.73 


11. 3(11. 5)^^ 


SSC2 


3.34x10-" (1%) 


2300.0(104.7) 


< 3.4 


3.0 


3.7x105 


0.65 


10.4(9.7)'^ 


SSCS 


5.77x10-15 (2.5%) 


787.0(53.4) 


6.8(2.5) 


7.0 


7.1x105 


0.92 


29.1(26.6)'! 


SSC4 


2.42x10-15 (3.9%) 


176.0(10.7) 


12.4(1.7) 


11.0 


1.1x10^ 


0.20 


32.9(26.9)° 


SSC5 


2.29x10-15 (4.5%) 


84(4.4) 


15.1(2.3) 


13.0 


2.9x10'^ 


0.84 


38.6(29.5)'' 


SSC6 


6.34x10-16 (7.7%) 


187(22.5) 


13.9(1.9) 


11.0 


2.6x105 


1.08 


20.6(21.49)"^ 



Note. — The first column shows the Ho fluxes and the EW(Ha) determined in the present work. 
The errors are shown inside the brackets. The ages in the fourth and fifth column are obtained from 
the EW(Ha) by R08 and this work, respectively. The uncertainties associated to our age estimates 
are ^ 1 Myr. The following two columns show the values of the masses and visual extinctions 
produced by the model fit, at given age, of the full SEDs of each cluster. In the last column, we 
listed the predicted radii of the Hii around the clusters. See text for details. 

''Estimated ages by R08 comparing the observed EW(Ha) with SB99 predictions. 

''Ages estimated in this work using the re-estimated EW(Ha) and Cloudy predictions for our models. 

"uh = 10* cm-3 

'^UH = 10^ cm-3 

"UH = 5 x 10^ cm-3 



TABLE 2 

Ages, masses, and extinctions of the 
clusters produced by the SED FIT. In this 
case age is treated as a free parameter. 



Ids Age [Myr] Mass [IO^Mq] A^, [mag] 



SSCl^ 


o n+0.0 
ci.U_2 Q 


4.7 


0.73 


SSC2'' 


s.otro 


3.7 


0.65 


SSC3 




5.1 


0.85 


SSC4 


Q n+2.0 
^•^-3.0 


8.8 


0.20 


SSCS 


y-U-5.0 


28.0 


0.81 


SSC6 




2.9 


0,93 



Note. — 

''The best fitting ages, masses, and extinctions are 
the same as the one derived in the first set of fits 
for the youngest SSCl and 2. See text for details. 



tron temperature of 20000 K (determined by Izotov et 
al. (1997) from optical emission lines of the galaxy). 

The values we estimated are listed in Table 1. We de- 
rived quite small sizes for the Hii regions of SSCl and 2, 
due to the high value of uh of the photoionized cloud sur- 
rounding the recently born clusters. The sizes increase 
for lower values of nn- This behavior is expected since 
the ionizing photons travel farther in less dense environ- 
ments. Nevertheless we found consistency between the 
assumptions made in the models and the aperture used 
by R08 to make the photometric measurement. 

5. DISCUSSION 

5.1. The origin of the IR excess in the SSCs 

R08 presented evidence of a flux excess in the red 
optical and NIR wavebands. The excess at ~ 0.8 /xm 
was attributed to the ERE phenomenon. The flux ex- 
cess at wavebands between 1.6 and 2.1 //ni was con- 
sidered to have two different origins. For tlie youngest 
clusters SSCl and SSC2, hot dust (~ 800 K) emission 



was advocated to explain the rise of the IR continuum. 
For the remaining SSCS, 4, 5 and 6, they showed that 
the infrared colors were consistent with red super giants 
(RSGs) which mainly contribute to the NIR light at ages 
above 7 Myr. 

We propose a new interpretation of the flux excess in 
the SBS 0335-052E SSCs. We have shown that the inclu- 
sion of the emission from photoionized gas in the models 
can well reproduced the colors (Figure 4) and the SED 
shape (Figure 5) of all the SSCs from UV to NIR. The 
ages found using the new estimated EW(IIq:) are very 
similar to the one previously determined from RGB and 
confirm the general picture of the star formation prop- 
agating from North to South (Thompson et al. 2006; 
R08). Moreover previous spectro-photometric studies of 
the galaxy from optical to infrared confirmed the pres- 
ence of dominant nebular contribution (Vanzi et al. 2000, 
Izotov et al. 2001). Vanzi et al. (2000) estimated that 
more than 45 % of flux in the NIR spectra of the SSCl-l-2 
and SSC4+5 regions are due to nebular continuum and 
emission. They also found from H/3/Br7 ratio a visual 
extinction of Ay « 0.73, which is in good agreement 
with the values we found from the SED flts. Finally, 
they interpreted the red H ^ K color observed globally 
in the galaxy as due to a hot dust component (^670 
K). However due to the very bright Br7 emission and 
low extinction observed across the two knots (SSCl-l-2 
and SSC4-I-5), they concluded that even if there may be 
a large quantity of dust in the galaxy, it is not located 
inside the SSCs. Due to the very young ages of SSCl 
and SSC2, the absence of any treatment of the nebu- 
lar emission, and the very high temperature of the dust 
advocated to explain the red excess in R08, make this ex- 
planation difficult to reconcile with the previously found 
observational properties. We showed that the integrated 
fluxes of these two clusters can easily be fltted from UV 
to NIR with self-consistent models including both stellar 
and nebular components. 

For the remaining clusters R08 showed that the colors 



8 



Adamo, A. et al. 



10 



10 



E 10 



-15 



-16 



10 



10 



10 



-19 - 



0.5 



SSC 1 1 

observed SED i 

n UVtoNIR tit (Stars + neb): 

UV-UBV tit (Stars only)- 




1.0 1.5 2.0 

A(/xm) 



.-1 4 [ 



10' 



-15 



observed SED 
n UVtoNIR tit (Stars + neb): 
g UV-UBV tit (Stars only)- 



E 1( 



-16. 



,-17 



10 




JL- 



1-15 



0.5 



1.0 1.5 

A(/im) 



2.0 




0.5 



SSC 5 

observed SED 

□ UVtoNIR tit (Stars + neb) 

UV-UBV fit (Stars only) 



1.0 1.5 



2.0 



SSC 6 

observed SED 
□ UVtoNIR tit (Stars + neb) 
UV-UBV tit (Stars only) 




2.0 



A(/im) 



Fig. 5. — SED analysis of the 6 clusters. The filled blues points are the observed integrated fluxes and associated error bars (R08). The 
black squares show the best fitted integrated fluxes produced by the our models. We plotted also the corresponding spectrum in black. 
In red we show the best flt produced by stellar continuum models. The red diamonds show the data points which have been included in 
the fit. The red crosses show the integrated fluxes excluded from the fit. In all the cases modeled spectra and integrated fiuxes have been 
corrected for the estimated extinction. See the text for more details. 



are compatible with the presence of RSGs, which are 
important in the NIR range at cluster ages around 7—10 
Myr. 

However, Papaderos et al. (2006) found for SSC3 evi- 
dence of Wolf-Rayet stars (WRs) in the optical spectrum. 
The presence of WRs is a very transient phenomenon and 
Starburst99 models predict that they last around 1 Myr 
at a cluster ages of around 6 Myr. This age is in agree- 
ment with our estimate for SSC3. It is likely that SSC3 
is still too young to show contribution due to RSGs. 



In the NIR spectrum of the knots SSC4-I-5, Vanzi et 
al. (2000) found no evidence of CO absorption typical 
of RSGs and important if the nebular emission does not 
contribute any longer (McCrady et al. 2003). We have 
to mention that the regions where SSC4 and SSG5 are 
located are very complex due to the proximity of very 
young star forming sources. It is possible that the spectra 
presented by Vanzi et al. (2000) are indeed dominated by 
these very young star forming regions. 

However we detect Ha emission at the location of 
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Fig. 6. — Evolution in time of the QLyc produced by the Star- 
burst99 model for a 10® cluster. Ovcrplotted the estimated 
QLyc of the six SSCs from the measured, extinction corrected 
f(Ha). The error bars show uncertainties of 25 % on the esti- 
mations. 

SSCS, SSC4, SSCS, and SSC6, which reveal that the 
emission from the photoionized gas is not gone yet. We 
have showed in Section 3 that the fraction of nebular 
emission contributing in the NIR filters (bottom panel 
of Figure 3) is still considerable and is strictly connected 
with the metallicity of the gas. The same models but 
at higher metallicity (i.e. solar metallicity) show that 
the nebular phase is over at around 6 Myr. Since SBS 
0335-052E is among the known lowest metallicity nearby 
galaxies, we expect to see such prolonged contribution. 

Both R08 and T09 analysed the Paa NICMOS data 
available for SBS 0335-052E, reaching discordant results 
regarding the Paa emission in the SSC3, SSC4, SSC5, 
and SSC6. T09 reported a none detection of Paa emis- 
sion in these clusters, showing a dearth of ionised gas. 
On the other hand, the values reported by R08 (listed 
in their Table 3) are close to the ones predicted by the 
Ha /Paa ratio considering our derived extinction. 

5.2. What is the cause of the I band excess? 

In the previous section, wc noticed that the six clusters 
have the observed flux in the I band (F791W) that sit 
above the best model predictions, with a discrepancy of 
0.2 — 0.5 mag. It is likely that the flux excess in this fll- 
ter is due to another mechanism than photoionized gas. 
The presence of a strong UV radiation field fuelled by 
the very young massive stars and large quantities of dust 
distributed in the galaxy and close to the SSCs can orig- 
inate the ERE phenomenon as already proposed by R08. 

An other possible explanation for this small systematic 
offset could be related to the low resolution power of the 
F791Wdata. The coarse pixel scale of the WFPC2/WFS 
chip (0.1") used to capture the F791W image and the 
small aperture radius used (0.15") means that the aper- 
ture correction is significantly more uncertain than for 
the remainder of the UV/optical data which was ob- 
tained with ACS. A significant error in the F791W aper- 
ture correction would lead to the SSCs lying systemati- 
cally above or below the best fit models. The accuracy 
of the encircled energy distribution determination for a 
0.15" radius aperture is no better than 0.07 magnitudes 
(Holtzman et al. 1995). While this is smaller than the 
observed excess (~ 0.3mag) it cannot be excluded that 



an uncertainty in the aperture correction contributes to 
the apparent offsets. The clusters are poorly sampled 
in the image and this could produce higher uncertainties 
associated with the derived aperture correction. 

6. CONCLUSIONS 

We tested our models wliic;li include stellar continua 
and a self-consistent treatment of the photoionized gas 
to verify the impact of the nebular component on the 
total flux of the young star clusters. Using new mea- 
surements of the EW(q;) we determined the ages of the 
6 bright SSCs in SBS 0335-052E. The ages were used to 
produced new fits to the whole SED of the clusters using 
our custom-designed evolutionary models. With the con- 
tribution to the integrated fluxes from the photoionized 
gas, we can fairly reproduce the SED shape in all the 
six clusters. This new proposed explanation is also sup- 
ported by previous numerous analysis of the starburst 
regions in the galaxy. 

We verified whether the assumptions made by our 
models are in agreement with the observed physical prop- 
erties of the clusters and in general, of the starburst en- 
vironment. We also tested the condition that the whole 
ionized region contributing to the total integrated fluxes 
is contained in the aperture size used to do photometry. 
We found that the size of the emitting Hii regions around 
the clusters are reasonably smaller than the photometric 
radii. 

We finally noticed a persistent flux excess at 0.8 fim, 
impossible to reconcile with our spectral evolutionary 
models. We consider, in agreement with R08, the pos- 
sibility that dust photolumincscence, typically observed 
in Hii regions around massive star forming complexes, 
can contribute to the integrated flux transmitted by the 
WFPC2/F79m filter. 

While nebular emission seems to be the origin of the 
red excess in the two dwarf galaxies SBS 0335-052E (this 
work) and NGC 4449 (Reines ct al. 2010, 2008a), this is 
not the case for the red excess observed in the young star 
clusters of Haro 11 (Adamo et al. 2010), which remains 
unresolved. The access to X-shooter data for the SBS 
0335-052E SSCl+2 and few of the clusters with strong 
red excess in Haro 11, would represent a fundamental 
source of information to disentangle between the different 
causes of the red excess, and quantify the real amount of 
flux produced by nebular continuum and emission. 
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